Wetting transitions in polydisperse fluids 
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The properties of the coexisting bulk gas and liquid phases of a polydisperse fluid depend not only 
on the prevailing temperature, but also on the overall parent density. As a result, a polydisperse fluid 
near a wall will exhibit density-driven wetting transitions inside the coexistence region. We propose 
a likely topology for the wetting phase diagram, which we test using Monte Carlo simulations of a 
model polydisperse fluid at an attractive wall, tracing the wetting line inside the cloud curve and 
identifying the relationship to prewetting. 
PACS numbers: 64.70Fx, 68.35.Rh 



Complex fluids such as colloidal dispersions, polymers 
and liquid crystals are inherently polydisperse, that is 
their constituent particles exhibit continuous variation 
in terms of some attribute such as their size, shape or 
charge. The presence of polydispersity in these sys- 
tems can profoundly influence their thermodynamic and 
processing properties [1-3]. Recent work has shown 
that the bulk phase behavior of polydisperse fluids can 
be surprisingly rich, exhibiting novel features such as 
polydispersity- induced phase transitions in polymers [4] , 
reentrant melting in colloids [5], and a sensitivity of co- 
existence properties to the presence of rare large parti- 
cles [6] . This richness is traceable to fractionation effects 
[3]: at coexistence the distribution of the polydisperse 
attribute varies from phase to phase. 

Notwithstanding the progress in elucidating the role of 
polydispersity in determining bulk phase behaviour, lit- 
tle is known regarding its effects on the wetting behavior 
of complex fluids. Besides being a matter of great funda- 
mental interest (see e.g. [7-9]) an understanding of this 
behaviour is crucial for technological applications such as 
spin coating and ink-jet printing of polymer or colloidal 
films for devices and displays [10]. 

This Letter broaches the subject of wetting in polydis- 
perse fluids. Starting from heuristic arguments based on 
the restricted case of a binary mixture, we propose likely 
wetting scenarios for polydisperse mixtures. These we 
test using detailed Monte Carlo simulations of a model 
fluid at an attractive wall. Our results reveal a novel form 
of wetting transition occurring inside the coexistence re- 
gion, which is driven by changes in density and/or tem- 
perature, rather than by temperature alone as in the case 
of monodisperse systems. We study the locus of wetting 
transitions in the phase diagram and how it depends on 
the strength of the wall-fluid attraction. For strongly at- 
tractive walls we further identify a prewetting transition 
within the stable vapor phase region. 

The state of a polydisperse system (or any of its 
phases) is described by a density distribution p{(j), with 
p{a)da the number density of particles in the range of 
a . . .a + da; for e.g. size polydispersity, a would be the 
particle diameter. In the most commonly encountered ex- 
perimental situation, the form of the overall or "parent" 



distribution, p*'(cr), is fixed by the synthesis of the fluid, 
and only its scale can vary depending on the proportion 
of the sample volume occupied by solvent. One can then 
write /o"(cr) = n°/°(cr) where /"(u) is the normalized 
parent shape function and n° — N/V the overall particle 
number density; as the latter is varied, p^{o') traces out a 
"dilution line" in density distribution space. The phase 
diagram is then spanned by n'^ and the temperature T. 

Phase separation in a polydisperse system occurs 
within a region of n*^ and T delineated by the cloud curve, 
which records the temperature T where phase separation 
first occurs for given n*^. Quite generally, one finds that at 
two phase coexistence the parent p° (a) will split into two 
"daughter" phase distributions p^((t) and p^{a) with dif- 
ferent shapes — this is the phenomenon of fractionation. 
The daughter distributions are related to the parent via a 
simple volumetric average: {i — £,)P^{'^)~^^p'^{'^) = P^{'^)^ 
where 1 — ^ and ^ are the fractional volumes of the two 
phases. In contrast to a monodisperse system where the 
densities of the coexisting phases are specified by the bin- 
odal and depend solely on temperature, full specification 
of the coexistence properties of a polydisperse system re- 
quires not only knowledge of the cloud curve, but also 
the dependence of ^, p^(cr) and /9^(cr) on n° and T [11]. 

To understand the effects of polydispersity on wetting 
transitions, we turn first to the simplest case of a bi- 
nary mixture with densities pi and p2 of two species 
of particles of different sizes; these are the analog of 
p{<j). Fig. 1(a) sketches an exemplary [12] isothermal 
cut through a bulk phase diagram, showing a region of 
liquid-vapor phase separation with tie-lines that shrink 
to zero at a critical point (CP). We assume that the 
particle-wall interaction strength (averaged over the two 
species) is sufficient that there is a finite region of coex- 
istence, bordering the CP, for which the liquid wets the 
wall (the wall is always wet at criticality because of the 
rapid decrease of the liquid- vapor surface tension [13]); 
a prewetting line extending into the one-phase region is 
also indicated. 

The dilution line constraint of a fixed shape for p^{a) 
reduces in the bidisperse case to a fixed ratio p\l p\- As 
n° — P1+P2 is increased from zero, the system follows the 
dilution line into the coexistence region. For a given point 
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FIG. 1: Schematic phase diagram of the binary mixture de- 
scribed in the text. Thick solid line: coexistence boundary; 
thin solid lines: tie-lines; grey areas: liquid (dash-dotted line) 
wets wall; dashed line: dilution lino; dotted line: prcwctting 
line, with CP (open circle). Full circle: bulk CP. (a) and 
(b): isothermal cuts at two different T as indicated in the full 
phase diagram (c) for fixed p\l p\. 



on the dilution line, the parent splits into two daughter 
phases located at the ends of the tie-line which intersects 
this point; owing to fractionation, the daughters lie o/f 
the dilution line. In the example of Fig. 1(a), one sees 
that as the system enters the coexistence region the liquid 
does not initially wet the wall, but as tl^ incrc^ascs fur- 
ther, the properties of the daughter phases change, and 
it starts to do so; crucially, this change occurs within the 
coexistence region. Fig. 1(b) shows a possible scenario at 
a higher temperature; here the phase coexistence region 
has shrunk and thus the dilution line now crosses the 
prewetting line in the one-phase region, indicating that a 
prewetting transition will occur on increasing n°. In this 
case, the liquid daughters always wet the wall, and so no 
further transitions occur inside the coexistence region. 
Generalizing from these two fixed temperatures to a full 
nP — T phase diagram, we infer the topology shown in 
Fig. 1(c); the wetting transitions inside the coexistence 
region will occur for a range of temperatures and so de- 
fine a wetting line which we expect to link the low and 
high density sides of the coexistence region. Provided 
this wcrtting line is first order where it intersects the gas 
branch, it will connect to a prewetting line, as shown. 

We have not, of course, explored all possible phase 
diagram structures here; oven for bidispcrsc systems at 
constant T, other possibilities could be envisaged. We are 
also not aware of any rigorous argument that indicates 
that the wetting line must have a negative slope as in our 
illustration. However, at least for weak wall interactions, 



where the wetting line is expected to enclose only a small 
region around the CP, its slope should follow that of the 
cloud curve there; in all the examples of which we are 
aware, this slope is negative, and so we may anticipate 
that the scenario we have outlined is quite general. 

To test the above ideas in a gemiincly polydispcrse set- 
ting, we have deployed Monte Carlo simulation to study 
a model fluid in which the polydisperse attribute a is the 
particle diameter. Particles interact via a Lennard-Jones 
potential of the form: 



(1) 



with dj = aiCTje/a"^, where e and a set our energy and 
length scale units respectively, while (Ty = {ui + o'j)/2 
and rij = |rj — rj|. The potential was truncated for 
Tij > 2.5aij and no tail corrections were applied. The 
diameters a are drawn from a (parental) Schulz distri- 
bution /°(cr) cx cr^ exp [{z + 1)<t/<t], with mean diameter 
(T. We chose z = 50, corresponding to a moderate degree 
of polydispersity: the standard deviation of particle sizes 
is 5 = l/^/z'+T « 14% of the mean. The distribution 
/°((j) was limited to within the range 0.5a < a < lAa. 

The bulk liquid-vapor phase behaviour of the model 
of eq. (1) has previously been obtained [6] within the 
grand canonical ensemble (GCE) using bespoke compu- 
tational techniques designed to ensure that the ensem- 
ble averaged density distribution always equals a desired 
parent form, by controlling an imposed chemical poten- 
tial distribution /i((T). Also implemented was a recently 
proposed strategy [11] which ensures that the two phases 
appear with equal probabilities in the GCE simulations. 
This results in finite-size corrections to the limiting bulk 
properties which are exponentially small in the system 
size - a feature which, it transpires, is crucial for accu- 
rately obtaining the wetting properties corresponding to 
a semi-infinite system. The study of ref. [6] yielded (in- 
ter alia) accurate estimates for the cloud curve and the 
form of /i(cr) at all points within the coexistence region. 
The CP was found to be located not at the temperature 
maximum of the cloud curve, but at higher density and 
lower temperature, as is also observed in experimental 
studies of polydisperse systems [14]. 

Equipped with knowledge of the bulk phase behavior 
of our model, we turn now to address its wetting prop- 
erties [13] by introducing a single attractive wall to our 
system. Our simulation box has dimensions 
with Lx = Ly = Iba and = 30a or Lz = 40a, and 
is periodic in the x and y directions. A fluid particle 
interacts with the attractive wall at 2: = via a potential 
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{ui + u)/2. The wall at 



with tiw = tw^cTi/a'^ and Uiw 
z = Lz'is, completely hard. 

In a canonical setting, the wall is "wet" at a given 
coexistence state point if the free energy is minimized 
when liquid is in contact with the wall, rather than gas. 
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In our GCE simulations, wc locate wetting points for a 
given wall strength by imposing the /i(cr) appropri- 
ate to some prescribed bulk values of and T; this 
corresponds to the conditions of a semi-infinite geome- 
try. We then record the ensemble averaged density profile 
n{z) = J dap{z,a) and the probability distribution 
of the fluctuating total density n = J da dz p{z,a); both 
of these are accumulated as histograms, one over z and 
the other over n. We repeat these measurements, for a 
range of bulk state points in which n° is increased isother- 
mally from the stable vapor region to values within the 
bulk coexistence region. Formally, wetting occurs within 
the GCE when the adsorption r{n^) — J^{n{z) — ny)dz 
diverges, with n„ the vapor density. For systems of finite 
size, the wetting transition is identifiable as the value of 
n° at which the instantaneous profile n{z) first becomes 
freely fluctuating. This leads to a linear decay of the 
n{z) along the .z-direction; it is also signalled by the dis- 
appearance of a "vapor peak" in the distribution of the 
overall density p(n) [15]. 

As an independent consistency check on this approach 
for determining wetting points, we have, for a broad se- 
lection of our state points, also implemented the distinct 
method described in ref. [16]. Here direct measurements 
(in geometries different to that considered above) of 
liquid- vapor, vapor- wall and liquid- wall surface tensions, 
together with Young's equation, j^y cos = ^vw~1lWi 
serve to estimate the contact angle 6. Wetting points 
follow as those combinations of rfi and T where 6 van- 
ishes. The requisite surface tensions were extracted from 
the form of p{n) (via multicanonical simulations [17]) for 
both a bulk (periodic) system and a system with two 
identical attractive walls (at 2; = and z = L^) as de- 
scribed in ref. [16]. In all cases, good agreement was 
found with the estimates of wetting points obtained from 
the observations of n(z) and p{n), described above, indi- 
cating that the effect of the hard wall si z = Lz on the 
freely fluctuating form of n{z) does not engender signifi- 
cant finite-size corrections to estimates of wetting points. 

Fig. 2 shows our measurements of the locus of wetting 
transitions in the n° — T plane, together with the cloud 
curve, for three choices of the wall strength ew Consid- 
ering first the weakest wall, = 1.04, we find that the 
line of wetting points traverses the coexistence region, in- 
tersecting the cloud curve discontinuously at its high and 
low density branches. The line has a negative gradient 
and encloses a portion of the coexistence region contain- 
ing the CP. For coexistence points inside this region the 
wall is wet, while for points outside, it is dry. Fig. 3 shows 
density profiles n{z) for two state points slightly on ei- 
ther side of the wetting line. When the wall is dry, the 
profile is gas like, except for a small enhancement at the 
wall. However, when the wall is wet, the liquid-gas inter- 
face fluctuates freely and the ensemble averaged profile is 
sloped because the maximum extent of the fiuctuations is 
limited by the hard wall a.t z = L^- For a system infinite 
in the z direction, the time averaged profile of the wet 
wall would become flat as a function of z. We have veri- 
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FIG. 2: Wetting phase diagram for three wall strengths. Pull 
line: cloud curve obtained in Ref. [6]. Symbols: measured 
points on the wetting line; dashed Unes guide the eye; small 
solid squares: state points referred to in Fig. 3. Inset: prewet- 
ting line observed for tm/t = 1.74. 



fied that the profile indeed becomes flatter as we increase 
Lz (inset Fig. 3), although no accompanying systematic 
changes in the locus of the wetting line were observed. 




FIG. 3: Profiles n{z) for the state points indicated in Fig. 2, 
see text. The inset shows the effect of increasing Lz on the 
profile of the freely fluctuating film. 

For the intermediate wall strength e^, = 1.21, the re- 
gion of coexistence for which the wall is wet appears to be 
larger than for the weak wall, although we were unable to 
follow the wetting line right across the coexistence region 
due to the computational burden of simulating at very 
low T and large nP . For both the weak and intermediate 
wall strengths, no sign of prewetting was discernible in 
the stable vapor phase region. However a clear signature 
of prewetting was found for the strongest wall strength, 
e-u, = 1.74, at low T in the stable vapor phase region near 
where the wetting line (of which we were able to deter- 
mine only a small segment) intersects the gas branch of 
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the cloud curve. The inset of Fig. 2 depicts the measured 
hne of prewetting transitions, obtained from the points at 
which the film thickness jumps from a thin to thick finite 
value, as shown in Fig. 4. We note that the prewetting 
line appears to intersect the cloud curve tangentially, as 
is mandated thermodynamically [18]. 

Although the prewetting line represents a true (sur- 
face) phase transition, we note that no fractionation ef- 
fects (such as splitting of the prewetting line into separate 
cloud and "shadow" curves - see Ref. [3]) are expected 
because neither of the coexisting phases are 'macroscopic' 
in the sense that their volume scales like the system vol- 
ume. We further note that the observation of prewetting 
is evidence of a first order wetting transition [13], while 
its absence at smaller wall strengths could indicate that 
critical wetting occurs here. Intriguingly, there seem no 
grounds to rule out, a priori, a scenario in which the 
wetting transition changes from first order to continuous 
wetting at some point along the wetting line, the two 
regimes being separated by a tricritical wetting point. 
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FIG. 4: Typical profiles n{z) for showing coexisting thin and 
thick films at the prewetting transition for Lz = 30a, see text. 



The new feature of the wetting phase diagram of Fig. 2, 
compared to its monodisperse counterpart, is the oc- 
currence of wetting transitions at constant T driven by 
changes in parent density. The wetting line delineates 
a subregion of the coexistence region for which the wall 
is wet, and which contains the critical point. We expect 
this subregion to shrink on,to the critical point as the wall 
strength is progressively weakened. We note that one can 
view the wetting line in the — T plane as the analog of 
the wetting temperature isotherm T = Ty, (a horizontal 
line) in the n — T phase diagram of a monodisperse sys- 
tem. The variation of daughter distributions with parent 
density in a polydisperse system causes this line to 'tilt'. 
The simulation results confirm our heuristic arguments 
to this effect based on a simple bidisperse example. 

Finally, with regard to the wider impact of our find- 
ings, we expect these to be relevant to the surface phase 
behavior of all complex fluids that exhibit significant 
polydispersity. More specifically, the possibility of den- 
sity driven wetting transitions could aid an understand- 
ing of the factors controlling the morphology of spin- 
coated or ink-jet printed thin films, which arises from the 
interplay of wetting and phase separation [10]. During 
these production processes, the solvent is progressively 
removed, thereby increasing the solute density, and pos- 
sibly leading to the type of wetting transitions we have 
described. 
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